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ABSTRACT 

Black hole masses are estimated for radio-loud quasars using several self-consistent 
scaling relationships based on emission-line widths and continuum luminosities. The 
emission lines used, H/3, Mgn A2798, and Civ A1549, have different dependencies 
on orientation as estimated by radio core dominance. We compare differences in the 
log of black hole masses estimated from different emission lines and show that they 
depend on radio core dominance in the sense that core-dominated, jet-on objects 
have systematically smaller H/3 and Mg II determined masses compared to those from 
Civ, while lobe-dominated edge-on objects have systematically larger H/3 and Mgn 
determined masses compared to those from C iv. The effect is consistent with the H/3 
line width, and to a lesser extent that of Mgn, being dependent upon orientation 
in the sense of a axisymmetric velocity field plus a projection effect. The size of the 
effect is nearly an order of magnitude in black hole mass going from one extreme 
orientation to the other. We find that radio spectral index is a good proxy for radio 
core dominance and repeating this analysis with radio spectral index yields similar 
results. Accounting for orientation could in principle significantly reduce the scatter 
in black hole mass scaling relationships, and we quantify and offer a correction for this 
effect cast in terms of radio core dominance and radio spectral index. 
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1 INTRODUCTION 

Black hole mass is an important fundamental quasar param- 
eter that can be measured directly via reverberation map- 
ping, whic h probes the motions of gas in the vicinity of the 
black hole (|Blandford fc McKedfl98a iPetersonll 19931 ). If the 
motion of gas in the broad line region (BLR) is dominated by 
the gravitational potential of the black hole, then the central 
mass can be estimated from the scaled virial product, 



Mbh = f 



AV 2 R 
G ' 



(1) 



where G is the gravitational constant, R is the distance 
to the broad line region gas that has velocity dispersion AV, 
and / is a scale factor of order unity that contains unknown 
information about the BLR geometry and orientation ef- 
fects. 

The factor / may vary between objects and possi- 
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bly even between emission lines in the same object, but 
determining / for individual objects requires an indepen- 
dent measure of the black hole mass. Instead, a measure 
of the statistical scale factor, < / >, is borrowed from the 
Mbh — o~* relationship for quiescent galaxies in order to cal- 
ibrate the zero-point of the reverberation mapping masses. 
One recent measure ment of the statist ical scale factor yields 
< / >= 5.25 ± 1.21 (|Woo et al.ll2010r i. 

Reverberation mapping then measures the central black 
hole mass by combining emission-line width {AV) and the 
size of the BLR (Rblr), inferred from the time delay in 
an emission line (Rblr = cr). At present, reverberation 
mapping black hole masses have been measu red for approx- 
imately 50 active galact i c nuclei (AGN) (e .g. JPeterson et all 
|2004 iBentz et al.ll2009l : iGrier et ajjlioia ) . 

Another important result has followed from reverber- 
ation mapping: the radius-luminosity (R — L) relationship. 
The time lag for emission lines to respond to var iation in the 
continuum is longer in more luminous objects (jKaspi et al.l 
12001 [2005L such that R blr oc L 1 ^ 2 when host galaxy con- 
tamination is removed (jBentz et al.l 12009 ). This result is 
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based on measurements of H/3, but li mited C iv-based re sults 
are consistent with this relationship (|Kaspi et alj|2007h . 

The R — L relationship anchors the single-epoch "mass 
scaling relationships," which provide an indirect estimate of 
black hole mass. The mass scaling relationships are widely 
used because they can be applied to any object with a single- 
epoch optical/ultra- violet (UV) spectrum measurement of 
a broad emission-line width and continuum luminosity, the 
proxy for the BLR radius. The velocity dispersion is usu- 
ally measured, for H/3, Mgn A2798, or Civ A1549, either 
from the full width at half maximum (FWHM) or the sec- 
ond moment of the line, oi ine , called the line dispersion. The 
two line measurements have different char acteristics and ad- 
vantages, thou gh au„e may be preferred (j Collin et al . 2006; 
iPetersonllioill ). We focus on the FWHM line measurement 
in order to facilitate comparisons with older studies and save 
an investigation of au ne for future work. 

Multiple studies have shown that the H/3 FWHM mea- 
surements used to estimate black hole mass depend on ori- 
entation. Radio core dominance (R), the ratio of the radio 
core flux density to the extended radio lobe flux density 
at 5 GHz rest-frame, is related to the angle between the 
radio axis and th e line of sight in the beaming model for ra- 
dio s o urces (e.g., Scheuer fc R eadhcad 1979; Qrr fc Browne! 



Il982l ; iGhisellini et al.lll993h . IWills fc Browne! (|l986l ) found 
that the FWHM of H/3 was significantly correlated with ra- 
dio core dominance such that there is an absence of very 
broad widths for objects with large core dominance. Their 
result is consistent with the motion of the H/3-emitting gas 
being primarily in a plan e perpendicular to the radio axis. 
iMcLure fc Dunlod (120011 ) test such a disc model for the BLR 
and find that it is a good fit to observations for H/3. 

The relationship between broad-line FWHM and core 
dominance found for H/3 has also been investigated for other 
broad lines, sometimes using radio spectral index as a proxy 
for radio core dominance. Radio spectr al index and ra- 
dio core dominance ar e well correlated (|Brothertonl 1 19961 ; 
iJarvis fc McLu^l2006T ) but it should be noted that radio 
spectral index is neither the same as radi o core dominance 
nor a perfect measure of orient ation (e.g., DiPompeo et al.l 
|2012| ). IJarvis fc McLurd (|2006h and lFine et all (|201ll ) find 



that the FWHM of Mgn h as an orientation dependence 
weaker than that of H/3 and lFine et al.l (|201ll ) find no ori- 
entation dependence for the FWHM of Civ, where radio 
spectral index is used to indicate orient ation. This result 
for Ci v is consistent with the result from I Vestergaard et al.l 
(2000) that FWHM Civ does not correlate with radio core 
dominance. 

Given the orientation dependence of FWHM for the 
broad emission lines used to estimate black hole mass, 
it has been suggested that black hole estimates based on 
single -epoch scaling relationships will depend on orie ntation 
(e.g.. iBaskin fc Laorl 120051 ; IJarvis fc McLurd l2006h . Black 
hole masses canno t be as accurately estimated without AV 
jAssef et al.ll2012l) . so we investigate the dependence of black 
hole mass estimates on radio c ore dominance for t he radio- 
loud (RL) subsample from the IShang et al.l |201lf ) spectral 
energy distribution (SED) atlas. 

In Section [2] we provide sample data and discuss spec- 
tral fitting and black hole mass estimates. Section confirms 
the correlations between FWHM and radio core dominance 
and presents the correlations between differences in the log 



of black hole mass estimates and radio core dominance. Here 
we also provide corrections to black hole mass based on the 
radio core dominance. To conclude Section the analysis is 
repeated for radio spectral index. The correlations between 
differences in the log of black hole mass estimated from dif- 
ferent emission lines and radio core dominance are discussed 
in Section[4]and Section[5]summarizes our investigation and 
corrections to black hole mass. 

Black hole masses and luminosities are calculated using 
a cosmology with Ho = 70 km s _1 Mpc -1 , J1a = 0.7, and 

n m = 0.3. 

Because we employ several values of "R" in this study, 
we explicitly state the definition of each parameter here to 
avoid confusion. Radio core dominance, R, is the ratio of 
the core flux density to the lobe flux density at 5 GHz rest 
frame. The ratio of the core luminosity at 5 GHz rest frame 
to the optical V- band luminosity, Ry, is an other measure of 
core dominance |Wills fc Brothertonll 19951 ). Radio loudness, 
R* , is the ratio of the total rad io luminosity at 5 GHz rest 
frame to the optical luminosity (jStocke et al.lfl992T ). The ra- 
dio spectral index, a (/„ ~ is a slope between two or 
more radio frequencies. 



2 SAMPLE, DATA, AND MEASUREMENTS 

For this investiga tion, we the R L subsamp l e of the 
Shang et al l (l201l|) SED atlas. See IWills et all j 19951 ) or 



Netzer et al.l (1 19951 ) for additional sample details. Many of 



the blazars originally included in this sample have been ex- 
cluded from this i nvestigation. 4C 11.69 , which was identi- 
fied as a blazar bv lRunnoe et all (|2012bl ) based on infrared 
variability, has also been excluded. Our final sample has 52 
objects. 

Originally assembled to study orientation, sample mem- 
bers were selected to have similar extended radio luminosity, 
which is thought to be isotropic. Thus, objects in the sam- 
ple are similar to each other but viewed at different angles 
indicated by their radio core luminosities. While this selec- 
tion was compromised by the original spherical aberration 
problem with Hubble Space Telescope, requiring selection of 
brighter targets, the extended radio luminosities are similar 
with a limited range versus log R (Figure[TJ left), and varia- 
tion in radio core luminosity is therefore largely responsible 
for variation in log R (right). 

We calculate radio core dominance, the ratio of the 
flat-spectrum core to the steep-spectrum extended lobes at 
5 GHz rest frame, for the sample. The radio data from the 
IShang et al l (|201lh SEDs are insufficient for this calculation, 
so we assembled radio core and total flux densities from high 
quality surveys and the literature. We applied a k-correction 
to the radio fluxes assuming slopes of a = and a — —0.7 
for the core and lobes, respectively. We also calculate the 
orientation indicator log Ry, the 5 GHz core flux density to 
the optical V-band flux density. 

We calculate radio spectral index for our sample in or- 
der to provide an orientation measure that may be easier to 
obtain and can be measured in some cases when log R can- 
not. We use the convention /„ ~ v a , where /„ is the radio 
flux, v is the frequency, and a is the radio spectral index. We 
measure the slope in two ways from the radio data in the 
IShang et al.l (|201ll ) SEDs. a a ii is the radio spectral index 
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Figure 1. Core and extended radio luminosities versus log R. The success of selecting objects with similar extended luminosities is 
apparent in the left panel, so variation in log R is due primarily to variation in the radio core luminosity seen in the right panel. 



measured by fitting a line to all of the available radio data 
in logarithmic space. ot\'_g is the radio spectral index mea- 
sured from two data points at observed-frame 1.4 GHz and 
4.9 GHz. We include this second measure in order to avoid 
the turnover of the radio spectrum around 1 GHz that is 
due to synchrotron self absorption and does not reflect the 
orientation of the object. We do note that we do not observe 
such a turnover in most of our objects; most are well fit with 
a single line using all of the available radio data from the 
SEDs. 

Table [T] lists radio properties for the sample. It is ar- 
ranged as follows: 

Column (1) gives the object name. 

Column (2) gives the redshift from lShang et all (|201 if ) . 
Column (3) gives the frequency at which fluxes were ob- 
served for the core and total radio components, respectively. 
In cases where the extended radio flux was observed instead 
of the total, the frequency of this observation is listed in 
Column (3). 

Columns (4) and (5) give the k-corrected, 5 GHz rest- 
frame flux densities in mjy. In most cases, the extended flux 
is calculated from measurements of the core and total flux, 
but in some cases the extended and core emission were the 
quantities directly listed in the reference. The cases where 
extended emission is actually observed rather than calcu- 
lated are indicated in Column 5. 

Column (6) gives the log of the k-corrected, rest-frame, 
radio core dominance, log R. 

Column (7) gives the r adio loudness parameter, 
log R*, where IShang et al.1 (|201ll ) calculates R* = 
/(5 GHz)//(4215 A) in the rest frame. 

Column (8) gives the reference for the radio data used 
to calculate radio core dominance. Radio spec tral index is 
alway s calculated from the radio portion of the lShang et al.1 
|201ll ) SEDs. 

Spectral fitting to measure emission-line widths and 



contin uum luminosities was carried out by iTang et all 
(2012) following the prescription of IShang et al.l |2007). 
I Tang et al.1 fit the opti cal/UV region of the SEDs us- 
ing SPECFIT |Krissl [l994). They model the regions around 
Civ A1549, Mgn A2798, and R/3 each with a power-law 
continuum and a line profile of two Gaussian components. 
A third Gaussian, that is not included in the total FWHM 
of the line profile, is added to the H/3 profile to fit the con- 
tribution to emission from the narrow-line region. The C IV 
line does n ot have a strong c ontribution from the narrow- 
line region (| Wills et all 1 19931 ) and, while Mgn does some- 
times have such a component, it is not apparent in our ob- 
jects. Optical and UV Fell emission is modeled with tem - 
plates (|Boroson fc Greenlll992j ; IVestergaard fe Wilkesll200ll ) 
based on the Seyfert 1 galaxy, I Zwl, that are allowed to 
vary in amplitude and velocity width. Continuum luminosi- 
ties based on flux m easurements from these fits are from 
iRunnoe et all (1201231 ). The resulting FWHM values for H/3, 
Civ, and Mgn, the bolometric luminosity, and other rele- 
vant data are given in Tabled 

ITang et all (|2012l ) estimate black hole masses for the 
IShang et alj atlas using singl e -epoc h scaling relationships 
from [Vestcrgaard & Peterson! 1120061) for H/3 and C iv and 
from IVestergaard fc OsmeJ (|2009h for Mgn. Multiple lines 
are used to statistically approach the true black hole mass 
and these scaling relationshi ps were chosen specifical ly be- 
cause they are self-consistent (Vestcrga ard et ai1l201ll ). thus 
enabling meaningful comparisons among them. The opti- 
cal/UV region of our SEDs facilitates consistent estimates 
of black hole mass from different broad emission lines be- 
cause they each have quasi-simultaneous data and coverage 
over a large enough wavelength range to include H/3, Mgn, 
and C iv. Black hole masses are listed in Table [5] 

Not all measurements are available for all objects. 
There are 38 (73%) objects with measurements of log R, 
FWHMh/3, and FWHM c ,v, 39 (75%) objects with mea- 
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Table 1. Radio data 



Ohirrt 


Redshift 


"core/ ^tot 
(GHz1/(GHz1 

I VT I ±Zj J / \ VT.LJ.Zj I 


f 

L COT& 

(mjy) 


(mjy) 


lop- R 




L V.l...lt_.± l... Ill_A_ 




Q 4.9 


3C 110 


0.7749 


4.85/ 4.85 


33.38 


33.63 


-0.25 


2.68 


8 


-0.99 


-0.65 


3C 175 


0.7693 


8.47/ 4.86 


32.34 


34.17 


-1.84 


2.99 


2.9 


-1.02 


-1.10 


3C 186 


1.0630 


1.70/ 1.70 


34.19 


34.12 


0.06 


3.33 


11 


-1.20 


-1.18 


3C 207 


0.6797 


4.86/ 1.40 


33.81 


34.20 


-0.39 


3.63 


2,12 


-0.67 


-0.65 


3C 215 


0.4108 


4.89/ 4.89 


31.84 


33.36 


-1.52 


3.37 


6 


-0.99 


-1.00 


3C 232 


0.5297 


8.55/ 8.55 


33.64 


33.00 


0.64 


2.87 


7 


-0.51 


-0.41 


3C 254 


0.7363 


5.00/ 5.00 


32.45 


34.17 


-1.72 


3.71 


11 


-1.06 


-1.25 


3C 263 


0.6464 


5.00/ 4.85 


33.24 


34.16 


-0.91 


3.00 


2,9 


-0.87 


-0.90 


3C 277.1 


0.3199 


1.70/ 1.70 


31.98 


33.41 


-1.43 


3.53 


11 


-0.76 


-0.87 


3C 281 


0.6017 


1.40/ 1.40 


32.48 


33.44 


-0.97 


3.23 


10 


-0.88 


-0.85 


3C 288.1 


0.9631 


5.00/ 5.00 


31.99 


34.16 


-2.17 


3.42 


11 


-0.96 


-1.07 


3C 334 


0.5553 


1.40/ 1.40 


33.18 


33.83 


-0.65 


3.11 


12 


-0.87 


-0.83 


3C 37 


0.6661 


5.00/ 5.00 


33.20 


33.63 


-0.43 


3.74 


10 


-0.85 


-0.79 


3C 446 


1.4040 


1.40/ 1.40 


35.54 


34.01 6 


1.53 


4.34 


1 


-0.32 


-0.11 


3C 47 


0.4250 


4.89/ 4.89 


32.53 


33.81 


-1.29 


3.82 


6 


-0.96 


-0.93 


4C 01.04 


0.2634 


5.00/ 5.00 


32.48 


33.10 


-0.62 


3.41 


13 


-0.54 


-0.48 


4C 06.69 


1.0002 


1.40/ 1.40 


34.87 


32. 77 6 


2.10 


3.32 


1 


0.16 


0.37 


4C 10.06 


0.4075 


4.99/ 4.99 


32.64 


33.30 


-0.66 


2.67 


18 


-0.79 


-0.96 


4C 12.40 


0.6836 


5.00/ 5.00 


31.93 


33.09 6 


-1.16 


3.03 


16 


-0.97 


-0.99 


4C 19.44 


0.7192 


1.40/ 1.40 


34.29 


33.86 


0.43 


3.42 


12 


-0.44 


-0.03 


4C 20.24 


1.1135 


1.40/ 1.40 


34.14 


34.64 


-0.50 


3.62 


12 


-0.52 


-0.39 


4C 22.26 


0.9760 


1.40/ 1.40 


32.93 


33.95 


-1.03 


3.26 


12 


-0.92 


-1.03 


4C 30.25 


1.0610 


1.40/ 1.40 


32.16 


33.90 


-1.74 


2.92 


12 


-0.91 


-1.07 


4C 31.63 


0.2952 


1.40/ 1.40 


33.52 


32.54 


0.98 


2.93 


1 


-0.15 


0.00 


4C 39.25 


0.6946 


15.00/ 15.00 


33.46 


35.68 6 


-2.22 


3.65 


3 


0.08 


0.75 


4C 40.24 


1.2520 


15.00/ 15.00 


34.60 


35.08 


-0.48 


3.94 


3 


-0.19 


-0.01 


4C 41.21 


0.6124 


1.40/ 1.40 


33.40 


33.88 


-0.48 


2.91 


12 


-0.79 


0.36 


4C 49.22 


0.3333 


0.40/ 0.40 


32.99 


33.04 


-0.05 


3.36 


3 


-0.64 


-0.60 


4C 55.17 


0.8990 


1.40/ 1.40 


34.77 


33.53 


1.24 


3.74 


12 


-0.32 


-0.44 


4C 58.29 


1.3740 


1.40/ 1.40 


32.70 


34.26 


-1.56 


2.66 


12 


-0.99 


-1.15 


4C 64.15 


1.3000 


5.00/ 5.00 


32.48 


34.27 


-1.79 


3.37 


11 


-0.67 


-0.77 


4C 73.18 


0.3027 


1.40/ 1.40 


33.86 


32.60 


1.26 


3.20 


1 


-0.14 


-0.06 


B2 0742+31 


0.4616 


1.40/ 1.40 


33.52 


33.16 6 


0.36 


2.77 


12 


-0.60 


-0.28 


B2 1351+31 


1.3260 


5.00/ 5.00 


33.22 


33.70 


-0.47 


2.95 


17 


-0.79 


-0.80 


B2 1555+33 


0.9420 


1.40/ 1.40 


33.25 


33.04 


0.21 


2.99 


12 


-0.52 


-0.55 


B2 1611+34 


1.3945 


15.00/ 15.00 


35.22 


35.21 


0.01 


3.77 


3 


0.00 


-0.34 


MC2 0042+101 


0.5870 


1.40/ 1.40 


32.15 


32.98 


-0.82 


2.92 


10 


-0.87 


-0.78 


MC2 1146+111 


0.8614 


5.00/ 5.00 


31.89 


33.20 


-1.31 


2.55 


16 


-0.86 


-0.91 


OS 562 


0.7506 


15.00/ 15.00 


34.41 


33.52 


0.89 


3.35 


3 


-0.08 


0.24 


PG 1100+772 


0.3114 


4.89/ 4.89 


32.23 


33.32 


-1.10 


2.65 


6 


-0.84 


-1.00 


PG 1103-006 


0.4234 


1.40/ 1.40 


32.88 


33.27 


-0.39 


2.94 


12 


-0.67 


-0.55 


PG 1226+023 


0.1576 


1.40/ 1.40 


34.32 


33.68 


0.64 


3.22 


12 


-0.19 


-0.21 


PG 1545+210 


0.2642 


1.40/ 1.40 


31.50 


33.25 


-1.75 


3.00 


12 


-0.74 


-0.77 


PG 1704+608 


0.3730 


1.40/ 1.40 


32.27 


33.65 


-1.39 


2.82 


12 


-0.76 


-0.83 


PG 2251+113 


0.3253 


5.00/ 5.00 


30.72 


33.22 


-2.50 


2.46 


15 


-0.78 


-0.83 


PKS 0112-017 


1.3743 


15.00/ 15.00 


34.37 


34.93 


-0.57 


3.45 


3 


0.13 


0.23 


PKS 0403-13 


0.5700 


1.40/ 1.40 


34.55 


32.85 


1.70 


3.73 


1 


—0.36 


—0.35 


PKS 0859-14 


1.3320 


4.90/ 4.90 


34.94 


33.54 6 


1.40 


3.43 


8 


-0.24 


-0.29 


PKS 1127-14 


1.1870 


1.66/ 1.66 


35.29 


33.34 


1.95 


3.88 


4 


0.04 


-0.21 


PKS 1656+053 


0.8890 


5.00/ 5.00 


33.97 


33.99 6 


-0.02 


3.10 


5 


-0.19 


-0.00 


PKS 2216-03 


0.8993 


1.40/ 1.40 


34.57 


33.65 


0.92 


3.23 


1 


-0.08 


0.13 


TEX 1156+213 


0.3480 


1.40/ 1.40 


31.24 


32.47 6 


-1.23 


2.38 


12 


-0.75 


-0.79 



a Frequencies are observed-frame. Fluxes, log R, and R* are k-corrected to 5 GHz and in the rest-frame. 
" Extended radio flux was observed ra ther than calculated from to tal and core radio flu xes. 
c Ref eren ces: (lllCooper et alj d2007h: (2llMullin et all ll2006ri: (3llKovalev etldl d2005ll: (4llStanghellini et al.l 
ll2005h; (51 IScott et alj J2004I) ; f6l lGilbert et alj ll2004ll; ( 71 iFev fc CharlotJ feOOdl ; (8) iReid et all Jl999]l; (9 ) 



iDennett-Thorpe et alTlll997r ) ; f 101 iHutchings et alj [|l996h; (11) iReid et all <ll995l1 ; ( 121 iBecker et al.l 
(13l| Lister et alj d 1994h ; (14) ISpencer et alj dl99 ill; ( 151 iKellermann et al.1 lll989l l; (161 iHutchings et alj 
(17 nRogora et al l 1119891 ; (181 iMilev fc Hartsuiikei] ill97sl v 

** Radio spectral index measured from all radio data available in lShang et al. 1 1|2011|1 SEDs. 
e Radio spectral index measured between observed-frame 1.4 and 4.9 GHz. 



1995); 



1988); 



© 2012 RAS, MNRAS 000, [T]- 



Orientation dependence of quasar black hole masses 5 



surements of log R, FWHM H/ 3, and FWHM Mg n, and 51 
(98%) objects with measurements of log R, FWHMciv, and 
FWHM Mg ii. 



3 ANALYSIS 

We investigate the dependence of the difference between the 
log of the black hole masses estimated with different scaling 
relationships on radio core dominance. We first demonstrate 
the relationship between FWHM and radio core dominance 
for H/3, Mgn, and C IV, and then proceed to characterize the 
dependence this creates in black hole masses derived using 
scaling relationships calibrated for each of these emission 
lines. We then repeat the analysis for radio spectral index. 

3.1 The orientation dependence of broad-line 
widths 

The scatter between black hole masses estimated from dif- 
ferent emission lines that results from varying radio core 
dominance is derived largely from the broad-line FWHM 
dependence on radio core dominance. Figure [2] shows log R 
versus FWHM for all three broad emission lines, where log R 
is on the y-axis for an easy comparison with previous results. 
We confirm that H/3 has the strongest log R dependence, fol- 
lowed by Mg II, and finally C IV which shows no dependence 
on radio core dominance. The given correlation coefficients 
are Spearman rank correlation coefficients and the proba- 
bility of finding these distributions of points by chance are 
listed in each panel. 

3.2 The orientation dependence of black hole 

masses 

The radio core dominance dependence of broad-line FWHM 
for H/3 and Mgn implies that black hole masses esti- 
mated from those line s may also have a simila r dependence 
( Bas kin fc Laorl |2005| ; IJarvis fc McLurel |2006| ). iTang et al.l 
( 20121 ) show that there is scatter between black hole masses 
estimated from different broad emission lines. We quantify 
the scatter around the one-to-one relationship by taking the 
standard deviation of the residuals in the y direction from 
the best fit to the parameters: around H/3 and C IV is 0.40 
dex, around Mgn and Civ is 0.194 dex, and around H/3 
and Mgn is 0.20 dex. This scatter is in part due to radio 
core dominance and can be reduced by applying a correction 
based on log R. 

3.2.1 Correlations 

In order to investigate how different black hole mass scal- 
ing relationships might depend on radio core dominance, 
we correlate the difference between the log of black hole 
masses estimated from different scaling relationships with 
radio core dominance. Table [3] gives the Spearman rank cor- 
relation coefficients and corresponding two-tailed probabili- 
ties. We normalize the masses rather than simply correlating 
masses derived from H/3, Mgn, and Civ individually with 
radio core dominance because we want to isolate the bias 
in black hole mass from orientation and avoid effects associ- 
ated with the sample having a range in intrinsic black hole 



Table 3. Correlations for black hole mass, log R, and a a u 



Param 1 




Param 2 


IN 




n a. 

p 




Alog(MBjf)(H/3 


-Mg n) 


log R 


39 


-0 


392 


0.014 


Alog(M Bff )(H/3 


-Civ) 


log-R 


38 


-0 


515 


0.001 


Alog(M SH )(Mg 


ii-Civ) 


log R 


51 


-0 


321 


0.022 


Lcore. 




\ogR 


52 





837 


< io- 6 






log R 


52 


-0 


294 


0.034 


Alog(M M )(H/3 


-Mg n) 


Lext 


39 





313 


0.052 


Alog(M SH )(H/3 


-Civ) 


Lext 


38 





228 


0.168 


Alog(M si/ )(Mg 


ii-Civ) 


Lext 


51 


-0 


154 


0.282 


Alog(M sff )(H/3 


-Mgn) 


Lcore. 


39 


-0 


259 


0.111 


Alog(M Bff )(H/3 


-Civ) 


Lcore 


38 


-0 


435 


0.006 


Alog(M SH )(Mg 


n-Civ) 


Lcore. 


51 


-0 


449 


0.001 


Alog(M SH )(H/3 


-Mgn) 


Oiall 


39 


-0 


426 


0.007 


Alog(M Bff )(H/3 


-Civ) 


Otall 


38 


-0 


495 


0.002 


Alog(M sff )(Mg 


n-Civ) 


Oiall 


51 


-0 


340 


0.015 


log/? 




Oiall 


52 





625 


< io- 6 



a p is the Spearman rank correlation coefficient. 

^ P is the given probability that the observed distribution of 

points would arise by chance. 

mass, which may mask orientation effects. The caveat to this 
approach is that if any systematic differences between the 
line width measurements depend on orientation, they will 
be folded into our results. 

Table [3] indicates that the log R dependence of H/3 and 
Mgn FWHM is propagated to the black hole masses de- 
rived from the H,9 and Mgn scaling relationships. Taking 
into account the number of objects available for each corre- 
lation and the probability indicated by the Spearman rank 
correlation test, the difference between the log of black hole 
masses estimated from H/3 and Civ, Mgn and Civ, and 
H/3 and Mgn is correlated with radio core dominance with 
a t-ratio of 3.5, 2.4, and 2.7, respectively. This correlation 
with radio core dominance causes a range in the difference 
between black hole masses estimated from different emis- 
sion lines that covers nearly an order of magnitude from 
the most jet-on (large values of log R) to the most edge-on 
(small values of log R) objects. In practice, the actual range 
will depend on the range of orientations present in any given 
sample and uncertainty in the orientation from log R. 

Core and extended radio luminosities are also listed 
in Table [3] in order to verify the source of the corre- 
lation between H/3 and Mgll-derived masses and radio 
core dominance. The differences Alog(M_BH)(H/3— C iv) and 
Alog(A/BH )(Mgll— C iv) are more significantly correlated 
with core luminosity than the extended luminosity, indicat- 
ing that the correlations between the mass difference and 
log R are due to variation in the core luminosity among 
objects that are intrinsically similar. 

3.2.2 Regression analysis 

We also performed a regression analysis in order to describe 
the trends with radio core dominance. The fits in this pa- 
per are m ade by minimizin g the chi-squared statistic us- 
ing mpfit (Markwardt 2009) which employs the Levenberg- 
Marquardt least squares method. The uncertainties in log R 
are considered to be negligible compared to those in black 
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Table 2. Black hole masses and luminosities 

Object log(L bo 0^ FWHMh/3 
(ergs s -1 ) (km s" 1 ) 


FWHM Mg „ 
(km s- 1 ) 


FWHMctv 
(km s- 1 ) 


log M SH (H/3) 
M 


log M BW (Mgn) 

M 


log M SH (Civ) 
M Q 


9P i in 


a ft on 

4b.yy 


1 o/i ten 


5645 


5700 


10.09 


9.47 


9.47 


oL 1 / 


A ft Ol 




9515 


6915 


10.60 


9.97 


9.58 


QP I Oft 

oL^ loO 


A ft ftft 
40.00 




5790 


6290 




9.35 


9.35 


oL zU f 


A ft OQ 
40. ZO 


00U0 


4140 


4935 


8.74 


8.92 


8.92 


op Ol c 

OL/ zlO 




( OU 


6260 


5605 


9.04 


8.92 


8.72 


op non 


A ft Qft 
40. oO 


4D00 


4620 


7145 


9.11 


9.16 


9.35 


QP OC/I 

oL zo4 


/)(! 17 


1 /inntc 
14Uyo 


6595 


5205 


9.89 


9.27 


8.96 


QP OftQ 

oL> zoo 


/Ift Q 1 

40. ol 


/i (i7n 
4y / U 


4265 


3310 


9.31 


9.22 


8.94 


op O r 7'7 1 


/i K fti 
40.01 


OOOO 


3380 


3215 


8.41 


8.34 


8.24 


QP OQ1 

oL zol 


A ft OQ 

40. zo 


'zoo c; 
/ yoo 


5505 


4865 


9.38 


9.10 


8.93 


OQQ 1 


/i ft fti 
40.01 


Oy/U 


4320 


4015 


9.63 


9.06 


8.93 


QP QQ/1 

OLj oo4 


/Ift /I A 
40.44 


0o40 


4840 


5745 


9.31 


9.16 


9.24 


QP Q7 


A ft 1 Q 

40. lo 


/i oqo 
4zoU 


4200 


3360 


8.62 


8.72 


8.39 


QP A A P. 

oL 440 


/i 9 m 
4 f .Ul 




3955 


3390 




9.21 


8.87 


QP / 1'7 

oL 4/ 


/I c; f)7 

4o. y f 


1 Anne 
14UU0 


7745 


5450 


9.64 


9.17 


8.82 


4L> U1.U4 


40.44 


none; 
yyUO 


9550 


6665 


9.25 


9.14 


8.54 


a p nft ftn 
40 uo.oy 


A 1 QO 

4 f .oU 


A m c; 
4U10 


4220 


5620 


9.33 


9.38 


9.61 


/if 1 1 n oft 
4L/ 1U.U0 


/Ift /1Q 

40. 4o 


4 ; 00 


4670 


3785 


9.00 


9.05 


8.81 


/ir^ 1 o /in 
4L/ lz.4U 


/I ft 1 1 
40.1 f 


OOOO 


4205 


5300 


8.58 


8.76 


8.92 


/IP 1 Q A A 

4^ iy.44 


/ifi 01 

40. y 1 


40 / 


3105 


2730 


9.23 


8.92 


8.75 


/i C" 1 on oa 
4L> zU.z4 


a ft no 

40. yz 




4070 


3525 




9.14 


8.93 


zip on oft 
40 ZZ.ZO 


40. oy 




4910 


5015 




9.07 


9.04 


4L/ oU.zo 


A ft 07 

40. z f 




4440 


3730 




8.87 


8.67 


/IP Q1 ftQ 

40 ol.Oo 


/IK fi1 
40.01 


ooyo 


4055 


4840 


8.78 


9.07 


9.21 


A C~* QCi oc; 
4L/ oy.zo 


/i ft m 
40. yi 


04UU 


5075 


4775 


9.51 


9.33 


9.21 


zip /in O/i 
4l_> 4U.Z4 


/IK fin 
40. OU 




3650 


4920 




8.90 


9.08 


4L/ 41. zl 


/l ft vc; 
40. ( O 


QA A CC 

o440 


3060 


3800 


8.80 


8.81 


9.03 


/if 1 /in oo 
4U 4y.zz 


/I c; no 

40. yy 


oyiu 


4415 


4535 


8.51 


8.69 


8.61 


zip r: c; 17 
4L/ OO.l ( 


/Ift /Ift 
40.40 




3475 


6420 




8.82 


9.25 


/if 1 ctq on 
4L^ Oo.zy 


A 1 OO 
4 f .ZU 




4830 


5745 




9.50 


9.64 


/I f 1 PA 1 K 
40 04.10 


/I ft 

40. / 




5125 


7245 




9.29 


9.52 


A f 1 Ti 1 Q 


/ift /i 
40. 4U 


oUyO 


3345 


3560 


8.66 


8.81 


8.73 


rSz U f 4z-f-ol 


/ift /ift 
40.40 


tuoyu 


6060 


4890 


9.85 


9.40 


9.04 


T30 1 or:i 1 01 

r5z lool+ol 


/I ft ftft 
40.00 




4205 


3690 




9.06 


8.89 


do 1 c;cc: 1 00 
r5z 1000-t-oo 


AP Ol 
40. Zl 




4650 


4240 




8.91 


8.75 


T30 1 ftl 1 I 1A 
OZ lDll + o4 


A 1 Oft 

4 f .U0 


4 i yo 


4010 


4625 


9.41 


9.22 


9.32 


IVIOZ UU4Z-|-1U1 


40. Oo 


oz / u 


6275 


4195 


9.21 


8.97 


8.43 


A/rr^o 1 1 /ift 1111 
IVLL/z 1140-hlll 


/l ft 00 

40. zy 


700 c 


5020 


3715 


9.46 


9.10 


8.67 


PC CftO 
Uo OOz 


A ft V/l 
40. / 4 


00UO 


3465 


3470 


8.86 


8.93 


8.92 


tjp 1 1 nn i V70 
11UU+ f / z 


/ift /ift 
40.40 


yoyu 


6090 


4775 


9.57 


9.26 


9.03 


op 1 1 HQ nnft 
rli HUo-UUO 


A ft QO 

40. oU 


0Z / U 


4780 


4515 


9.00 


9.01 


8.85 


pp 1 ooft 1 noQ 
l^O IzzOH-Uzo 


A ft Oft 

40. yo 


O4U0 


2935 


4530 


8.96 


8.92 


9.27 


"DP 1 1 010 


/I HQ 

40. yo 


OftoO 


5230 


4560 


9.12 


8.94 


8.71 


PG 1704+608 


46.49 


10465 


8405 


4015 


9.81 


9.61 


8.87 


PG 2251+113 


46.35 


4060 


4110 


4805 


8.89 


8.96 


8.94 


PKS 0112-017 


46.88 




4130 


5030 




9.21 


9.36 


PKS 0403-13 


46.36 


3735 


3535 


3325 


8.82 


8.78 


8.60 


PKS 0859-14 


47.22 


4615 


3775 


4520 


9.46 


9.29 


9.40 


PKS 1127-14 


47.15 




3725 


3695 




9.18 


9.12 


PKS 1656+053 


47.06 


3510 


3980 




9.13 


9.22 




PKS 2216-03 


46.95 


4415 


3555 


3600 


9.35 


9.14 


9.04 


TEX 1156+213 


46.01 


7740 


5175 


3880 


9.19 


8.91 


8.56 



a From lRunnoe et al] l|2012ah . all other columns from lTang et al] ||2012| V 



hole mass because the fluxes are well known (although the 
uncertainty in 9 is large) unless otherwise noted. 

The results of the regression analysis also suggest that 
the H/3 black hole mass scaling relationship has the strongest 
dependance on radio core dominance, though Mgn also has 
a dependence, albeit weaker. Figures [3] [4] and [5] show three 



combinations of difference in log of black hole mass ver- 
sus log R with the regression line over-plotted. The differ- 
ence Alog(M bh)(R/3— Civ) shows the steepest slope versus 
log R, indicating that H/3 has the strongest bias with radio 
core dominance. The difference Alog(A/s_f/)(Mgn— Civ) has 
a slightly shallower slope, indicating that Mgll has a weaker 
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Figure 2. Radio core dominance versus FWHM for H/3, Mgll, and Civ. The effect seen in IWills fe Browne! jl986h . where large values 
of FWHM have small vales of log R but small values of FWHM have a range of log R, is reproduced for H/3 and, to some extent, Mg II. 
This effect is not statistically significant for C IV. Spearman rank correlation coefficients (p) are given with the probability (P) of finding 
this distribution of points by chance listed underneath. 




Figure 3. Log M B h (H/3)-log M sff (Mgn) versus log R for 39 
objects where measurements are available for all three quantities. 
The solid red line shows the linear regression with a nonzero in- 
tercept. The dependence on log R is weaker; the relationship is 
significant at nearly the 3 a level. 



Figure 4. Log M B ir(H/3)-log M flfr (Civ) versus log R for 38 
objects where measurements are available for all three quantities. 
The solid red line shows the linear regression with a nonzero in- 
tercept. The relationship is significant at greater than the 3 a 
level. 



dependence on log R than does H/3. The H/3 and Mgn- 
derived masses both depend on radio core dominance, but 
the negative slope of Alog(AfsH)(H/3— Mgll) versus log R 
supports the conclusion that H/3 has the stronger depen- 
dence. 

The significance of this effect can be increased by con- 
sidering the one-tailed probabilities for the correlation coef- 
ficients. Because we expect H/3 scaling relationships to un- 
derestimate black hole masses compared to C IV and Mg II 
scaling relationships for jet-on objects, we can halve the 



probabilities given in Table [3] According to the one-tailed 
probabilities, the difference between the log of black hole 
mass versus log R is correlated with a probability of finding 
the observed distribution of points by chance of: < 0.7% for 
H/3 and Mgll < 0.05% for H/3 and Civ and < 1.1% for 
Mgll and C IV. 
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Figure 5. Log MaH(Mgn)— log Mbh(Civ) versus log R for 51 
objects where measurements are available for all three quantities. 
The solid red line shows the linear regression with a nonzero in- 
tercept. The relationship is significant at greater than the 2 a 
level. 



3.2.3 Luminosity considerations 

In order to verify that the dependence of difference in the 
log of black hole masses from different scaling relationships 
on radio core dominance is not the result of a luminosity ef- 
fect which enters via the R — L relationship, we consider 
our sample in the luminosity-log R plane. We expect to 
see a slight correlation because, although the selection for 
similar extended luminosities means that we are looking at 
similar luminosity objects in principle, orientation effects 
are l ikely to introduce a trend into the observed luminosity 
(e.g. iNemmen fc Brothertonl l^OlO) because accretion discs 
are not isotropic emitters. Figure [6] shows an increase in 
AL A with radio core dominance for our objects. More face- 
on objects (larger log R values) have higher apparent lu- 
minosities, which will result in higher calculated black hole 
masses for these objects because luminosity goes to the one 
half power in the mass scaling relationships. We observe the 
INemmen fc Brothertonl (|2010D effect and the dependence of 
difference in black hole mass on radio core dominance exists 
in spite of it. 

In principle, continuum luminosity and FWHM must 
both be corrected for inclination and an empirical correc- 
tion to black hole mass will incorporate both dependancies, 
though the primary effect is derived from the FWHM depen- 
dence, which enters as the square in black hole mass scaling 
relationships while luminosity is only to the one half power. 

In practice, we find that the empirical mass correction 
does not account for the luminosity dependence because we 
fit the correction based on the differences in the log of the 
black hole mass, so that the luminosity also enters as the 
difference, AL(5100A - 1450A) and AL(3000A - 1450A). 



Figure [6] shows that the optical and UV luminosities have 
nearly the same dependence on log R, so these differences 
are flat with radio core dominance. In order to remove 
the radio core dominance dependence from the luminosity 
for individual single-epoch mass estimates, we therefore 
provide corrections in Equations 12141 



At 1450 A: 



AL A 
AL A 



At 3000 A: 



al; 



ALa 
AL A 



(0.292 ±0.157) log R, 
(0.097 ±0.101) log R, 



(0.162 ±0.135) log R, 
(0.240 ± 0.104) log R, 



for log R ^ 
for log R < 



for log R ^ 
for log R < 



At 5100 A: 



al; = 



ALa 
ALa 



(0.218 ± 0.169) log R, for log R ^ 
(0.136 ±0.111) log R, forlogi?<0 



(2) 



(3) 



(4) 



ALa is the uncorrected monochromatic luminosity and 
AL; is the luminosity corrected for log R. The broken fit 
around log 7? = is naturally motivated by the shape of the 
data in Figure [B] 

The ratio of radio core luminosity at 5 GHz to a 
monochromatic optical luminosity, log Ry (different from 
radio loundess parameter, log R*) has been suggested to be 
a better indication of orientatio n than the radio core dom i- 
nance, l og R, with some caveats dWills fc Brothertonlll995T l. 
As with lVestergaard et al.l (|2000h . our result persists, but is 
slightly less significant, when performing the analysis with 
log Ry instead of log R. 



3.3 Correcting black hole masses for orientation 

The correlations presented demonstrate a likely orientation 
dependence that is propagated into black hole masses cal- 
culated from H/3 and Mg n FWHM measurements. The fact 
that no significant radio core dominance dependence is ob- 
served for the FWHM of C IV allows a correction for orienta- 
tion to H/3 and Mgll-derived masses. For H/3, we used the fit 
of Alog(A#gir)(H/3— C iv) versus log R to zero the difference 
between H/3 and Civ-derived black hole masses, normalizing 
at log R — 0. We note that other normalizations are possi- 
ble, and an edge-on normalization might be preferable but 
is not possible here as a true edge-on view is likely obscured 
by a dusty torus. With the orientation dependence removed 
from the difference, it can be combined with the orientation 
unbiased C IV black hole mass estimates to give an unbiased 
estimate of H/3 black hole mass. This process was repeated 
for Mgn. Equations [S] and [B] can be used to estimate black 
hole masses unbiased by orientation from H/3 or Mgn. 



lop 



log 



M M (H/3)' 



M, 



© 



Mbh (Mg ii)' 



log 



M BH (Kj3) 



© 



± (0.173 ±0.051) log R 



Al 



© 



= log 



M SH (Mgii) 



M 

+ (0.066 ± 0.028) log R 



(•») 



(6) 
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Figure 6. AL A at 1450, 3000, and 5100 A versus log R. The correlations between the difference in the log of black hole mass and log R 
occur in spite of the correlations between \L\ and log R. The Spearman rank correlation coefficient and corresponding probability of 
finding this distribution of points by chance are given. 



Mbh{B-P)' is the log of the H/3-based black hole mass 
with the log R dependence removed, and Mbh(H-/3) is the 
log of the orientation-biased, H/3-based black hole mass esti- 
mated from the scaling relationship. The same conventions 
apply for Mgn. The uncertainties given here are the \-a 
errors in the fit coefficients and are not independent. 

Applying this correction reduces the scatter between 
black hole masses estimated from different broad emission 
lines. The scatter around the one-to-one relationship is re- 
duced from 0.40 dex to 0.35 dex for Rf3 and C iv, 0.194 dex 
to 0.188 dex for Mgn and Civ, and 0.20 to 0.18 for H/3 
and Mgn. Figures [7] and [8] show the reduction in scatter 
that results from applying the corrections for log R for H/3 
and Mgn. Original masses are shown in the left panel and 
corrected masses are on the right. 

The caveat for these corrections is that they are linear 
fits to data that display a linear relationship in log space 
with large scatter. The errors are large due to the large scat- 
ter present in the data, and fitting with a different algorithm 
or prescription will yield a slightly different correction. 

3.4 Radio spectral index 

We repeat the analysis for radio spectral index, an orienta- 
tion indicator that is more easily measured than radio core 
dominance because the object does not need to be resolved 
in the radio image. We present only the results for radio 
spectral index measured from all the available radio data. 
The results for radio spectral index measured between 1.4 
and 4.9 GHz are similar, but the correlations are slightly 
weaker. 

Radio spectral index is a good proxy for radio core dom- 
inance, the two are correlated with a probability of less than 
10 -6 for the correlation to arise by chance. Figure [§] shows 
radio spectral index measured from all the available radio 
data points versus radio core dominance. We find a best- 
fitting line of: 

a a ii = (0.260 ± 0.048) log R - (0.433 ± 0.048) (7) 
The orientation dependence of FWHM of H/3 and Mg II 




log(R) 

Figure 9. Radio spectral index measured for all available radio 
data points versus radio core dominance. a a u is a good proxy for 
log R, which is traditionally used as an orientation indicator. The 
biggest outlier in this correlation has a spectral shape in the radio 
that is not well fit by a power law. 



and lack of a dependence for C IV that is indicated by a 
radio core dominance dependence is also present for radio 
spectral index. Figure \TU\ shows this dependence on radio 
spectral index for all three lines. 

In order to investigate how different black hole mass 
scaling relationships might depend on radio spectral index, 
we correlate it with the difference between the log of black 
hole masses estimated from different scaling relationships. 
The Spearman rank correlation coefficients and correspond- 
ing two-tailed probabilities are given in Table [3] 
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Figure 7. H/3 versus Civ-based black hole masses (left) and corrected H/3' versus Civ-based black hole masses (right). The dashed line 
is the one-to-one relationship. 
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Figure 8. Mgn versus Civ-based black hole masses (left) and corrected Mgn' versus Civ-based black hole masses (right). The dashed 
line is the one-to-one relationship. Note that, while the reduction in scatter is modest, the points now fall more evenly around the 
one-to-one line. 



Taking into account the number of objects available for 
each correlation and the probability indicated by the Spear- 
man rank correlation test, the difference between the log of 
black hole masses estimated from H/3 and Civ, Mgn and 
C IV, and H/3 and Mg II is correlated with radio spectral in- 
dex with a t-ratio of 3.3, 2.5, and 2.85, respectively. As with 
radio core dominance, the correlation with radio spectral 
index causes a range in the difference between black hole 
masses estimated from different emission lines that covers 
about an order of magnitude from the most jet-on (large 
values of a aH ) to the most edge-on (small values of a aH ) 



objects. Similarly, in practice, the actual range will depend 
on the range of orientations present in any given sample and 
uncertainty in the orientation from ot a ii- 

We also repeated the regression analysis in order to de- 
scribe the trends with radio spectral index. The uncertain- 
ties in a a ii are considered to be negligible compared to those 
in black hole mass. The trends are very similar to those for 
radio core dominance. Figures 1111 1121 and [13] show three 
combinations of difference in log of black hole mass versus 
O-aii with the regression line over-plotted. 

The behavior of the 1450, 3000, and 5100 A monochro- 
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Figure 10. Radio spectral index versus FWHM for H/3, Mgn, and Civ. The effect seen with radio core dominance is reproduced for 
radio spectral index. Spearman rank correlation coefficients are given with the probability of finding this distribution of points by chance 
listed underneath. 
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Figure 11. Log M-bh (H/3)— log Msjj(MgIl) versus radio spec- 
tral index for 39 objects where measurements are available for 
all three quantities. The solid red line shows the linear regres- 
sion with a nonzero intercept. The relationship is significant at 
approximately the 3 <r level. 



Figure 12. Log M^h (H/3) — log Mbh (Civ) versus radio spectral 
index for 38 objects where measurements are available for all three 
quantities. The solid red line shows the linear regression with a 
nonzero intercept. The relationship is significant at greater than 
the 3 a level. 



matic luminosities with radio spectral index is similar as 
for radio core dominance. The differences in luminosity are 
not correlated with a a ii but the individual luminosities are, 
so we provide a correction to these that should be applied 
before calculating black hole masses from the single-epoch 
scaling relationships. There is no evidence in Figure [14] that 
a broken fit is necessary here, so we fit only one line. The 
data have a lot of irregular scatter and we find that as- 
suming the errors in radio spectral index are negligible and 
weighting all the monochromatic luminosities similarly gives 



the best fit. We normalize the correction at a a U = —0.432, 
which corresponds to log R — in Equation [7] 
At 1450 A: 

AL' A = AL A - (0.716 ± 0.202) a a ii - (0.310 ±0.211) (8) 
At 3000 A: 

AL A = ALa - (0.678 ± 0.201) a a ii - (0.293 ± 0.210) (9) 
At 5100 A: 

AL' A = AL A - (0.675 ± 0.222) a a ii - (0.292 ± 0.229) (10) 
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Figure 14. AL^ at 1450, 3000, and 5100 A versus a a ;;. The correlations between black hole mass difference and a a u occur in spite of the 
correlations between XL\ and a jj. The Spearman rank correlation coefficient and corresponding probability of finding this distribution 
of points by chance are given. 
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Figure 13. Log Ms/j(Mgll)— log TsAbh{Civ) versus radio spec- 
tral index for 51 objects where measurements are available for 
all three quantities. The solid red line shows the linear regression 
with a nonzero intercept. The relationship is significant at nearly 
the 3 a level. 



corrections for black hole mass in terms of radio spectral 
index: 



+ (0.580 ± 0.170) a all 
(11) 



M BH (wy a ~ 




~Mbh(KP) 




= log 


Mq 



M BH (Mgu)' a 



(0.251 ±0.180) 



Al 







log 



M SH (Mgn) 



M. 







(0.247 ± 0. 



(0.107 ±0.106) 



Afg.fi (H/3)q is the log of the H/3-based black hole mass 
with the a a u dependence removed, and Msif (H/3) is the 
log of the orientation-biased, H/3-based black hole mass es- 
timated from the scaling relationship. The same conventions 
apply for Mgn. The uncertainties given here are the 1-a er- 
rors in the fit coefficients and are not independent. 

Applying this correction reduces the scatter between 
black hole masses estimated from different broad emission 
lines. The scatter around the one-to-one relationship is re- 
duced from 0.40 dex to 0.35 dex for H/3 and Civ, 0.194 dex 
to 0.188 dex for Mgn and Civ, and 0.20 to 0.18 for H/3 and 
Mgn. To the number of decimal places shown, the reduction 
in scatter is exactly the same as for log 7?. Figures [TS] and 
[16] show the reduction in scatter that results from applying 
the corrections for a a u for H/3 and Mgn. Original masses 
are shown in the left panel and corrected masses are on the 
right. 



ALa is the uncorrected monochromatic luminosity and 
XL'x is the luminosity corrected for a a u. 

The regression and correlation analysis with radio spec- 
tral index indicate that a likely orientation dependence is 
propagated into black hole masses calculated from H,9 and 
Mg II FWHM measurements and we again provide a correc- 
tion for this effect. We use the fits of Alog(M B _tf )(H/3-C iv) 
and Alog(A/s£r)(Mgll— Civ) to zero the a a u dependence, 
normalizing to a a u = —0.432. Using the same methodol- 
ogy as for radio core dominance, we arrive at the following 



4 DISCUSSION 

We correlated differences in the log of black hole mass esti- 
mates with the orientation indicators radio core dominance 
and radio spectral index. We found that the dependence 
of FWHM of H/3 and Mg II A2798 on orientation indicators 
propagates into black hole masses estimated from scaling 
relationships based on these emission lines, and derived cor- 
rections to H/3 and Mgn black hole masses for this effect. 
The orientation indicator dependence is strongest for the 
black hole masses based on measurements of H/3 FWHM, 



© 2012 RAS, MNRAS 000, [TV?? 



Orientation dependence of quasar black hole masses 13 




9 9.5 
log(M BH ) (C IV) 




9.5 

log(M BH ) (C IV) 



Figure 15. H/3 versus Civ-based black hole masses (left) and radio spectral index corrected versus Civ-based black hole masses 

(right). The dashed line is the one-to-one relationship. 
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Figure 16. Mgll versus Civ-based black hole masses (left) and radio spectral index corrected Mgn^ versus Civ-based black hole masses 
(right). The dashed line is the one-to-one relationship. Note that, while the reduction in scatter is modest, the points now fall more 
evenly around the one-to-one line. 



although differences associated with Mgn-based measure- 
ments also show some dependence. 

The range in the difference between black hole masses 
estimated from different emission lines covers nearly an or- 
der of magnitude from edge-on to jet-on, although it will 
depend on the orientation biases in any given sample. Note 
that this sample was originally selected to have a full range 
in log R and will likely have a different distribution of ori- 
entations than a flux limited or randomly selected sample. 

Because of the nature of the broad-lined objects in 
our sample, "edge-on" will likely be closer to 45° than 90° 



l|Barthe]||l98"gh . If quasars could be observed at a true edge- 
on inclination, we might expect to measure an even broader 
FWHM and estimate a correspondingly larger black hole 
mass. The scale factor / derived from such a mass might be 
closer to unity than the larger statistical values that are cur- 
rently available, implying that the current value is inflated 
to compensate for orientation and projection effects. 

Radio core dominance and radio spectral index are in- 
dicators of inclination angle and, while the measurement un- 
certainties on log R and a a ii are relatively low, converting 
to 9 will introduce additional scatter. The exact amount is 
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model dependent, but at least in the case of lGhisellini et all 
l| 19931 ). the scatter is substantial. The implication for this 
analysis is that the correction to H/3 and Mgn black hole 
masses were derived from a fitting procedure that assumed 
negligible uncertainty in log R compared to the uncertainty 
in black hole mass. A physical correction for inclination an- 
gle would require another fitting procedure that assumes 
uncertainty in both parameters. We provide only the cor- 
rections for log R and a a u as they are model independent 
and more useful in practice. 

Given the size and significance of the R-dependence of 
black hole mass scaling relationships, it would be advisable 
to recalibrate them with the radio core dominance depen- 
dence taken into account. The H/3 scaling relationship is 
based on reverberation mapping of the H/3 line, so this en- 
tails doing reverberation mapping on RL objects and also 
obtaining UV spectra and assuming that the BLR is the 
same in RL and radio-quiet objects. Ideally this would be 
done for a large sample of RL quasars with known orien- 
tations, although few of these have been done. These new 
orientation-independent single-epoch mass scaling relation- 
ships would be useful for the ~ 10% of quasars that are RL. 
They must also have high-quality radio maps for measuring 
log R which becomes more difficult with increasing redshift 
for a number of reasons. This could also be done with ra- 
dio spectral index which is easier to measure and can be 
measured for some objects in which log R cannot. 

While the C IV scaling relationship does not itself de- 
pend on orientation, its original calibration has scatter from 
orientation effects and should be corrected. C IV and Mg n 
scaling relationships are based on reverberation mapping 
of HB (jVestergaard fc Petersonll2006l ; IVestergaard fc Osmerl 
2009) that shows a radio core dominance dependence. Re- 
calibrating the C IV and Mg II scaling relationships with this 
effect taken into account should reduce their scatter. 

When black hole mass scaling relationships are re- 
derived, they should also be cast in terms of the line dis- 
persion, (June- The parameter an ne is the "best" line width 
measure for deriving black ho le mass in the sense that it 
gives the least biased result (|Petersonl 1201 ll ) . There is a 
larger range in FWHM than ou ne and, since the measure 
of line width is squared in scaling relationships, the extreme 
values possible for FWHM will introduce a larger range in 
black hole mass. If the black hole mass scaling relationships 
are re-cast in a consistent way for line dispersion as well as 
FWHM, then the (j; ine -based scaling relationships can be 
explored to see if they are more robust than FWHM against 
radio core dominance and radio spectral index. 

It has been argued that C IV may not be an appro- 
priate choice for estimating black hole mass from a single- 
epoch spectrum because there may be a strong wind com- 
ponent and the motions of the C IV gas may not be dom- 



inated by grav i ty (e . g.. iDunlod 12004: iBachev et al.l 12004 



IShemmer etail 12004 iBaskin fc Laorl 120051 ). This~is likely 
the case in many narrow-line Seyfert 1 galaxies that ex- 
hibit the characteristic triangular profile indicative of an 
outflowing wind and the black hole mass scaling relation- 
ships do not apply to these objects (none of the objects in 
this sample have such a line profile for Civ.) Quasars do 
not share these types of asymmetric profiles, but their pro- 
files are often blueshifted r elative to the systemic redshift 
of the object (| Wilked H984T ) . The objects with the largest 



C iv blueshifts exhibi t broader lines ( Brotherton et aklll994l ; 
iRichards et all 120021 : iHo et al.l 120121 ). but the correspond- 
ing increase in black hole mass is small, well within 
the uncertainty in the mass estimate. Despite these is- 
sues, it is not clear that C iv is any less trustworthy 
of a black hole mass indicator than H/3 or Mgn, which 
also have problems l|Vestergaard et all l201ll : lAssef et al.l 
l201ll ). R everberation mapping on Civ shows a virial rela- 
tions h ip dPeterson fc Wandelll 19991 . l2000l : IVestergaard et all 
l201ll : iPetersonl I26T1I T" and there are a number of advan- 
tages to using the Civ line (IVestergaard fc Peterson! 120061 ; 



I Warner et al.ll2003l ; IVestergaardl2002f r Most importantly for 
this study, C iv is not biased by orientation as measured by 
radio core dominance or radio spectral index as are other 
lines. The lack of a significant radio core dominance depen- 
dence for C iv FWHM may be due to the fact that there is 
a com ponent to the line that does not reverberate l|Dennevl 
|2012| ). if the reverberating component could be isolated one 
might expect it to behave more like the H/3 line. 

The correction for orientation, as indicated by radio 
core dominance and radio spectral index, is only one of a 
suite of corrections that will r educe the scatter in black 
hole mass scaling relationships. lAssef et al.l (|201ll ) perform 
a study similar to this one on a sample of 12 high-redshift, 
gravitationally lensed quasars. They specifically compare 
black hole masses based on Civ to those from H/3 and Ha 
and find that the residuals between the C iv-based masses 
and the H/3 or Ha-based masses correlate with the ra- 
tio of the UV to optical continuum luminosities. A cor- 
rection based on this color dependence reduces the scat- 
ter in the residuals by almost a factor of two. Additional 
sources of scatter in the black hole mass scaling relation- 
ships include the rol e of radiation pressu re in determining 
the line profile (e.g., Marconi et alj I2OO8I ) . Eddington ratio 



llOnken fc Kollmeier 



ability l|Woo fc Park 



2008), uncertainties from source vari- 
201ll ). and effects of low S/ N data and 



—31 

measurement uncertainties (|Dennev et al.ll2~01ll ). While our 
result is based only on RL quasars, it may be true of RQ 
quasars as well (as shown bv ljarvis fc McLurell2006l ). If ori 



entation indicators, like that of BorosonT i 20 1 if ) . "are devel- 
oped for RQ quasars, this effect can be investigated and 
possibly used to further reduce the scatter in scaling rela- 
tionships. 



5 CONCLUSIONS 



Using the RL subsample of the lShang et all l|201ll ) SED at- 
las, we have analyzed the dependence on radio core domi- 
nance and radio spectral index for differences in the log of 
black hole masses estimated from self-consistent scaling re- 
lationships and quasi-simultaneous optical-UV data for H/3, 
Mgn, and C iv. We correlated the difference between the log 
of black hole masses estimated from these lines with orien- 
tation indicators and found that H/3 and Mg 11-based masses 
are biased with respect to radio core dominance and radio 
spectral index as compared to masses estimated from C iv. 
The dependence is such that H/3 and Mgn will underesti- 
mate black hole mass compared to C iv for jet-on objects, 
the effect is stronger for H/3, and the range in the difference 
between black hole masses estimated from different emission 
lines covers nearly an order of magnitude in black hole mass 
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from our most jet-on to edge-on objects. We used the C IV 
FWHM-based mass estimates, unbiased with respect to ra- 
dio core dominance and radio spectral index, to derive the 
following corrections for black hole masses based on H/3 and 
Mg II FWHM in terms of log R and a all : 



that affect mass estimates and searching for new correla- 
tions, we can reduce the scatter in the scaling relationships 
even further. 



lop 



lop 



Mbh(H/3)' 



Mi 



lop 







Msjf(Mgn) 



M BH (R/3) 



M, 



© 



M, 



© 



+ (0.173 ±0.051) log R 
Mbh (Mg n) 



lop 



M 







+ (0.066 ± 0.028) log R 



(13) 



(14) 
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(0.251 ± 0.180) 

M B H (Mgn) 



log 



M, 





(0.107 ±0.106) 



(0.580 ±0.170) a aU 
(15) 

- (0.247 ± 0.088) a an 
(16) 



When the log R corrections are applied, scatter around 
the one-to-one relationship is reduced from 0.40 to 0.35 dex 
between H/3 and C iv-based black hole masses, 0.194 to 0.188 
dex for Mg II and C IV, and 0.20 to 0.18 dex for H/3 and Mg II. 
The scatter reduction is the same when using corrections 
based on radio spectral index. 

We investigated the role of luminosity in the corre- 
lation between differences in the log of black hole mass 
and orientation indicators and found that differences in 
luminosity do not depend on radio core dominance or radio 
spectral index and are therefore not being corrected in our 
empirical corrections. The individual luminosities at 1450, 
3000, and 5100 A do depend on the orientation indicators 
so we provide separate corrections that are to be applied to 
the luminosity before calculating black hole mass from the 
single-epoch scaling relationships: 



At 1450 A: 

r XL X - (0.292 ±0.157) log R, for log R >0 
I AL A - (0.097 ± 0.101) log R, for log R < 



ALA 



XL'x = AL A - (0.716 ±0.202) a aU - (0.310 ±0.211) (18) 
At 3000 A: 

AL , _ JaLa — (0.162 ± 0.135) log 7?, for logi? ^ 
A \ AL A - (0.240 ± 0.104) log R, for logi? < 

XL'x = AL A - (0.678 ± 0.201) a aU - (0.293 ± 0.210) (20) 
At 5100 A: 

AL , = f ALa - (0.218 ± 0.169) log R, for logi? ^ 
A \ ALa - (0.136 ±0.111) log R, forlogi?<0 

XL'x = AL A - (0.675 ± 0.222) a aU - (0.292 ± 0.229) (22) 

This type of correction is only one of a suite of many 
corrections that will tighten the black hole mass scaling rela- 
tionships. By tracking down additional physical parameters 
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